A PRIORI ESTIMATES FOR FLUID INTERFACE PROBLEMS 
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Abstract. We consider the regularity of an interface between two incompressible and inviscid 
fluids flows in the presence of surface tension. We obtain local in time estimates on the interface in 
H^^'^^ and the velocity fields in H^'' . These estimates are obtained using geometric considerations 
which show that the Kelvin-Helmholtz instabilities are a consequence of a curvature calculation. 



1. Introduction 

In this manuscript we consider the interface problem between two incompressible and inviscid 
fluids that occupy domains and 0^ in R", n > 2, at time t. We assume that M" = Qf U U St 
where St = dQf, and let p± : -^M, v± : — M", and the constant p± > denote the pressure, 
the velocity vector field, and the density respectively. On the interface St, we let N±{t,x), x (z St 
denote the unit outward normal of (thus iV+ + = 0), H{t,x) € {T^St)'^ denote the mean 
curvature vector, and k± = H ■ N±. We also assume that there is surface tension on the interface 
given by the mean curvature. Thus the free boundary problem for the Euler equation that we 
consider here is given by 



(E) 



I pivt + \/^v) = -Vp, X G \ 5t 

|v-z; = 0, xGM^xS'j, 

The boundary conditions for the interface evolution and the pressure are 



{dt + v±-V is tangent to 115'*^ 
p+{t,x) - p^{t,x) = K+{t,x), X £ St. 

where we introduced the notation v = v+Iq^ + V-ln_ : M" \ St ^ M", etc. 

The boundary conditions (BC) are a consequence of assuming that 1) the interface velocity is 
given by the normal component of the velocity = v± ■ N±, and that 2) the surface tension on 
the interface St is given by the mean curvature of the surface. A weak formulation for the Euler 
flow with this form of surface tension is given by 

f {pv)t + Vypv = -Vp + H{t, x)5{St), x€W 

where 5 is the Dirac mass distribution. This weak formulation implies the boundary condition for 
the pressure stated in (BC). 

Here we consider the problem where Q'^ is compact and derive a priori estimates, local in time, to 
the problem (E, BC) that prove bounds on v{t, •) G Fi'=(M" \ St) and St G hI''+^ for |A: > f + 1. 
The assumption that $7"*" is compact is not necessary. In fact the same estimates hold if we assume 
that St is either periodic or asymptotically flat. A more interesting observation is that our proof 
works verbatim for the case where there are several fluids occupying regions with interfaces SI- 
Thus we can treat more general setting than the existing literature. 
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The interface problem between two fluids has been studied extensively in the math and physics 
literature. In the absence of surface tension it is well known that the interface problem between 
two inviscid and incompressible fluids is ill-posed due to the Kelvin-Helmholtz instability, and it 
is argued on physical basis that the surface tension is a regularizing force that should make the 
problem well posed. In |BHL93j . Beal, Hou, and Lowengrub demonstrated that the surface tension 
makes the linearized problem well-posed. For the full nonlinear problem rigorous results have been 
obtained for irrotational velocities. In this case the problem can be completely reduced to the 
interface evolution with nonlocal operators. For this problem, Iguchi, Tanaka, and Tani |ITT97j 
proved the local well-posedness in 2 dimensions with initial interface almost flat and initial velocity 
almost zero. For the general irrotational problem, Ambrose |AM03j and, more recently, Ambrose 
and Masmoudi |AM06j proved the local well-posedness in 2 and 3 dimensions, respectively. 

We should note that for irrotational flows without surface tension the interface problem is given 
by he Birkhoff-Rott differential-integral equation. Several results were obtained in this case such 
as those obtained by Sulem, Sulem, Bardos, and Frisch |SSBF8lj . and Wu |Wu06j . We also note 
that without surface tension one can consider the Euler equation with discontinuous velocity fields. 
Although the interface problem is ill-posed, due to the Kelvin-Helmhotz instability, weak solutions 
to the Euler equation which may include such discontinuities have been considered by DiPerna and 
Maida |DM87) . Delort |De91j . and others. There is also a rich literature of numerical studies of the 
interface problem, see for example |HLS97j and references therein. 

A related problem to the interface problem is the water wave problem where there is only one 
fluid. In this case The Rayleigh- Taylor instability, instead of the Kelvin-Helmholtz instability, may 
occur. Such problems have been extensively studied and there is a vast literature on this subject, 
see for example, cWiTWj . |Wu99j . [Lin5] . 

Our approach in obtaining energy estimate for the interface problem is similar to the water waves 
problem treated in [SZOGj in that it is geometric in nature. It is based on the well known fact that 
these free boundary problems have a variational formulation on a subspace of volume preserving 
homeomorphisms. We use this variational approach to determine the terms that should be included 
in the energy. Of course these terms are identifled as being the highest order terms of the linearized 
problem to (E, BC). It is worth noting that from our analysis of the operators involved in the 
linearized problem, the Kelvin-Helmholtz instability appears naturally as a consequence of the 
negative semi- definiteness of the leading part of the unbounded curvature operator of the infinite 
dimensional manifold of admissible Lagrangian coordinate maps. The surface tension, created by 
the potential energy of the surface area, generates a higher order positive operator that makes the 
linear problem well-posed and help to establish the energy estimates. The well-posedness of the 
full problem well be addressed in a forthcoming article. 

Our paper is organized as follows. In section 2 we explain how to determine the pressure from the 
velocity. In section 3 we give a variational formulation of the problem as a constrained variational 
problem for volume preserving maps. We use this formulation to motivate our definition of energy. 
In section 4 we prove that our energy controls the Sobolev norm of the velocity and the mean 
curvature and derive bounds on the energy. Some of the details in the geometric calculations 
are omitted since they are given in details in SZ06 and are available as notes on the web at 
[http : / / WW . mat h . gatech . edu^zengch /notes/ notes 1 . p df . 

Notation All notations will be defined as they are introduced. In addition a list of symbols will be 
given at the end of the paper for a quick reference. The regularity of the domains is characterized 
by the local regularity of St as graphs. In general, an m-dimensional manifold Ai C M" is said to 
be of class or H^, s > y, if, locally in linear frames, M can be represented by graphs of or 
mappings, respectively. 

As in SZoej A^^ denote the inverse Laplacian with zero Dirichlet data, Tl± denote the harmonic 
extension of functions defined on St into 0,^, and 7V± denote the Dirichlet to Neuman operators in 
the domain fi^ . Given two fluids in fi^ with constant densities p± we denote hy Af = ^ A/"+ + ^7V_ , 
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the operator J\f^^ acts on function with mean zero and its range are also functions with mean zero. 
For any quantity q defined on M" \ St we write q = + 9-lf7_ where q± = (?1q±. 

2. Determining the pressure 

In this section we explain how to express the pressure in terms of the velocity in this setting 
which is less clear than the free boundary problem of water wave in vacuum where the boundary 
conditions of p± are obvious. 

To determine the boundary value of p± we take the dot product of Euler's equation (E) with N± 

-N± ■ Vp± = p±nt±{v± ■ N±) - p±v± ■ T>t±N± 
and using the fact that + =0, we obtain 

—Vn+P+ + —Vn^P- = v+ ■ Dt+iV+ + V- ■ Dt_7V_ - V„t_„t v^. 
P+ P- + - 

Substituting the formula for Y)t±N±, which has been calculated in |SZ06j . 

(2.1) Dt^A^i = -i{Dv±nN±))'^ 
we have 

— Vjv+P+ + — V7v„P- = Il+{vl,vl) + n_(t;I, vl) - 2V„T_„T v^, 
P+ P- + " 

where II-i- is the second fundamental form of St associated to A'^-i-, which satisfy n_|_ + n_ = 0. 

Since p± = 'H±{p±\st) + A^^Ap± in ilf, we have 

—M+P+ + —M-P- = — ^V7v+A7^Ap+ - — VAr_Al^Ap_ 
P+ P- P+ P- 

+ n+(r;J, vl) + n_(r;]:, vl) - 2V^t„^t on Sf 

The boundary condition p+ — P- = n+ on St stated in (BC) implies that on St 

11 1 

P± = AA^H N^K^ Vat, A:r^Ap+ VAr„Al^Ap_) 

Pt P+ P- 

+ Ii+{vl,vl) + n_(t;I,^I) - 2V,T_„T vi). 
Finally, since V • f = in M" \ St, we have from (E) 

(2.2) -Ap = p\/ ■ (V^v) = ptiiDv)"^, xGR^xSt. 
Therefore, 

P±\st = AA-i(-— A/VKT+ViV+ATHr(L't;)2 + V N^AzhiiDv)'^), 

(2.3) Pt 

+ U+ivl,vl)+n.{vI,vI)-2V,T_,T vi). 

One can verify that the quantity 7V^^ acts on in the above has zero mean on St and thus p is well 
defined by and (1^ . 

3. Lagrangian formulation and the energy 

This section is intended to explain the intuition behind the energy. It illustrates how to isolate 
the leading order nonlinear terms £/ and defined in ()3.22|) and (|3.25|) respectively. 

In his 1966 seminal paper |Ar66j . V. Arnold pointed out that the Euler equation for an incom- 
pressible inviscid fiuid can be viewed as the geodesic equation on the group of volume preserving 
diffeomorphisms. This point of view has been adopted and developed by several authors such as 
D. G. Ebin and G. Marsden |EM70j . A. Shnirelman |Sh85j . and Y. Brenier |Br99j . to mention a 
few, in their work on Euler's equations on fixed domains. It is this point of view that we adopted 
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to explain the motivation for our definition of energy for the water waves problem |SZ06j , and it is 
this same point of view that forms our starting point to determine the appropriate energy for the 
interface problem . 

3.1. Lagrangian formulation of the problem. Conservation of energy can be obtained from 
multiplying the Euler's equation (E) by v, integrating on \ St, and using (BC) to obtain the 
conserved energy Eq: 

(3.1) E^ = E^{Suv)= [ ^dx+ [ dS^ [ ^dx + S{St), 

JR" ^ JSt Jflt ^ 

where S{-) denotes the surface area. 

Let u±{t,y), y £ 17^, be the Lagrangian coordinate map solving 

(3.2) ^ = y(^t,x), x(0) = y, 

then we have v = ut o u~^, and for any vector field w on x G M" \ St, D^tD = (w o u)t o u^^ . 
Therefore in Lagrangian coordinates the Euler's equation takes the form 

(3.3) putt = -(Vp) o n, u(0) = idco, 

where the pressure p is given by ()2.2|) and (|2.3|) . 

Since v{t, •) is divergence free in M" \ St, then u±{t, ■) are volume preserving. Moreover, while 
ii+(t,-)|5p = ii_(t,-)|5(, may not hold, it is clear that ii+(t, 5o) = U-(t,So). Thus the Lagrangian 
coordinates maps satisfy: 

1) $± : f]± ^ ^±i^^) a volume preserving homeomorphism. 

2) s = d<^±{n^) = $(ao±) 

Define 

r = {$ = + ^-1^-; ^± satisfy 1 and 2 above}. 

As a manifold, the tangent space of F is given by divergence free vector fields with matching normal 
component in Eulerian coordinates: 

r$r = {u) : M" \ 50 ^ I V • It; = and u-:^ + w^\^(^So) = 0, where w = {w o ^>"^)}. 

Here as in we are following the convention that for any vector field X : <I>(r2o) ^ I^" its 

description in Lagrangian coordinates is given hy X = X o ^. 

Writing S{^) = /^(^q) ^5' for the surface area of ^{Sq), the energy Eq in Lagrangian coordinates 
can be written as: 

(3.4) Eo = Eo{u,Ut) = l [ p\ut\'^dy + Siu), iu,ut)eTr 

where the volume preserving property of u is used. This conservation of energy suggests: 1) TL be 
endowed with the L'^{pdy) metric*; and 2) the free boundary problem of the Euler's equation has 
a Lagrangian action 

/(n) = / / f^^^dydt - [ S{u)dt, u{t, ■) e T. 

Let denote the covariant derivative associated with the metric on T, then a critical path u{t,-) 
of / satisfies 

(3.5) ^tut + S'iu) = 0. 

In order to verify that the Lagrangian coordinate map u{t, •) satisfying (E) and (BC) is indeed 
a critical path of /, it is convenient to calculate ^ and S' by viewing L as a submanifold of the 
Hilbert space L^{W \So,pdy,W). 

'Including the density in the volume element pdy introduces a factor of - in front of the physical pressure. 
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(T$r)^ and orthogonal decomposition of vector fields. For any vector field X defined 
on <1>(M" \ 5o), Hodge decomposition suggests that we decompose X into X = w — Vip, with 
i> = i>+tn+ + i'-tn., so that w = w o <l> e T$r and Vip o <^ e (r^r)-^. For any Y G T^T, the 
orthogonality / Wip ■ Ypdy = implies 

p+tp+Y+ + p-i^-Y^dS = 0. 

Therefore, ip must satisfy p+ip+ = P-i>~ — on <^(5'o). This suggests that, for any $ G F, 

(r$F)^ = {-{S/iP) o $ I p+ijj^ = on $(So)}. 

To prove this claim, we only need to find such a ip given X. From X = w — Vip and w\_ + = 0, 
we have 

(3.6) X^+X^ = -Vn+^P+ - Vn^iP- = -M^p^ - Vtv+A^^AV' - Vjv_ Al^A^^. 

Since V • w = 0, we obtain 



(3.7) 



-Atp = V -X 

_^±\hso) = ^^"^ = -^-^-1 {X^ + X^- V^+A^^V • X - Viv_ AI^V • X) . 
It is easy to verify that w = X — Vip satisfies w = w o ^ ^T^T. 

Computing S!t and //$. Given a path u{t,-) G F and v = ut- Let St = u{t,So). Suppose 
tD(t, •) G r„(j)F, then the covariant derivative ^tiD and the second fundamental form II^(£^{w,v) 
satisfy 

wt = &tw + IIu{t){w,v), ^tw G T„(t)F, II^t){w,v) G (r„(i)F)-^. 

Let V = Ut o = V o and w = w o he the Eulerian coordinates description of ut and w, 
then for X = Dju> there exists p^^^ = +Pw v^n- '■ ^ ""(^j Sq) :— > M determined by (|3.7() 

such that 

(3.8) P+Pt,w = P-Pv,w on St, II{iv, v) = -(Vp^,^) o n G (r„(i)F)^. 
The Eulerian coordinates description of the covariant derivative is given by 

(3.9) &tw = i^tw) o n"^ = Btw + Vp^,^,. 

The terms involving X = Dtw in equation 1)3. 7() are expressed as follows. From wj: + if:': = on 
u(t, Sq) and identity (|2.1j) we have 

{'Dt+w+) ■ N+ + (Dt.w.) ■ iV_ = Dt^u;^ + T>t_w^ + V^tv+ ■ N+ + V^tv^ ■ iV_ 

=V,T_,Twi + V^T_Tt;i - n+(^T^u;T) _ U.ivl,wl). 

And since V • D^?/; = tr{DvDw) then p^^^^ is given by 

-Aptt,^^ = tv{DvDw) 



(3.11) 



A more useful way to express the boundary value p^ ^ is as follows. From the divergence decom- 
position formula 

(3.12) = V ■v± = V ■v1 + K±v^ + N±-Vn±v± on St, 
where V is the covariant derivative on St, implies 

(3.13) V^^v± ■ N± = V^TV± ■ iV± - K±w^v^ - V ■ {wivD + V^tu;^. 
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Thus, we have 

- VN+^+hT{DvDw) - VN_^Z^tT{DvDw)}. 
Moreover, for any smooth function / defined on St, we have from the Divergence theorem, 

/ -f\/N±^±^t-r{Dv±Dw±) dS = - Vfn± ■ V A^hT{Dv±Dw±) + fn±tT{Dv±Dw±)dx 
Jst Jnf 

Again, by the Divergence Theorem the first term integrates to zero and the second term can be 
written as 



-dS 



(3.15) 



/ -fVN^A^hriDv±Dw±)dS = [ -fV^^v± ■ N± + w^VM^ ■ v±c 
JSt JSt 

- I D'^fn±{v±,w±)dx. 
Jnf 

Thus, using the decomposition Vf'H± = V^/ + {J^±f )N± and letting / = —M~^g, we obtain 
(3.16) / gp^^dS= [ -w^v^{M++N'-)M-^gdS + [ D'^{n±{M-^g)){v,w)dx. 

J St ' J St JR"\St 

Computing S"(<I>). By the variation of surface area formula, for any id € T$r we have 

< S'{^),w >L^Rr^^So,pdy)= / i^+wi dS = K^W^dS. 

We need to find the unique representation S"(<I>) in T$r of the above functional. 

Lemma 3.1. For any smooth function /o : ^(5*0) — > M, let f± = ifc H±J\f^ ^-^t fo '^"■^ / 
+ J-Iq- , then we have V/ G T$r and for any id E T$r, 



fow^dS = w ■ Vfpdx. 

*{So) JR"\$(5o) 

The verification of the lemma is straightforward. Therefore, we have 

(3.17) y($) = Vp« where = ^^7^±A/'~W^K±. 

P+P- 

Splitting of the pressure. From (|2.3() . 1)3. 7() and (|3.17|) . it is clear that p{Pv,v+Pk) = P- Therefore, 
we obtain the well known equivalence between equation ()3.5() for critical paths of I and the Euler's 
equation (E) with the free boundary condition (BC). The Euler' equation can also be written as 

(3.18) Btv + Vp^,„ + Vp« = 0. 

Notice that the pressure splits into two terms, the first p^j^^ is the Lagrange multiplier, and the 
second p^ is due to surface tension. These two terms will be treated differently in the energy 
estimates. 

3.2. Linearization. In order to analyze the free boundary problems of the Euler's equation, it is 
natural to start with the linearization. The Lagrangian formulation provides a convenient frame 
work for this purpose. From (|3.5|) . the linearized equation is 

(3.19) ^^w + £{ut,w)ut + ^'^S{u){w) =0, iZ;(t, •) G T„(i_.)r, 

where ^ is the curvature tensor of the infinite dimensional manifold F. Below we calculate ^ and 
^'^S{u), which is a linear operator on T^F. Since these operators are self-adjoint, we will compute 
their quadratic forms. 

Computing ^'^S{u). The formula for ^^S{u) was given in |SZ06j 
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J St 

+ wi i-^stwi - w^iinp + v^TK±) ds, 



for any w S T„r, where V is the Riemannian connection and is the Beltrami-Lapalacian 
operator on St- Of course ^'^S{u){w, w) is independent of the choice of the + or — sign. 

Needless to say that this is a very comphcated expression for ^'^S{u) {€],€}). We will single out 
its leading order part. Since the value of ^'^S{u) does not depend on the choice of + or — sign, 
we compute with the + sign and assume that St is a sufficiently smooth hypersurface. From the 
Divergence Theorem, 

\^'^S{u){w,w)- [ \V^w^\'^dS\ 
Jst 

<\ i^+C^N+pt,,w + V^tu;+ • + V^tw^) dS\ + C\w\l2^s^y 

To estimate the integral on the right side, we use the splitting of Vjy+pt] w 

Vtv+P^,^ = J^+{pt,,w\st) + V7v+A;iAp+ = —Af+p^ - VN+A+hr{Dwf. 

P+ 

and identities (|l'{.15)) and (l^-i.!!-}!) to obtain 

\^^S{u){w,w) - [ \V^wi\'dS\ < I / p'^^^—M+K+dS\ + Ci\w\l,^s^^ + \w\l,^^„^st))- 
J St J St P+ 

Finally, from 1)3. 16() . we have 

(3.20) \^^S{u){w,w) - f |VV|' ^-^l < C{\w\l,^st) + HlHR-~.St))- 

J St 

Using Lemma l3. 11 we can define a self-adjoint positive semi- definite operator £/{u) on T^T as the 
leading order part of ^'^S{u). In the Eulerian coordinates, £/{u) takes the form 

(3.21) ^{u){w) = S7Utn++'^f-tn-, where f± = -^n±M-^M^{-As,)wi. 
Clearly £/ (u) is self- adjoint and satisfies 

(3.22) £/{iv,iv)= [ \V^w^\'^dS. 

Jst 

which is like a third order differential operator on R" \ St- From (|3.23|) . we can write 

(3.23) S!'^S{u) = (n) -|- at most 1st order diff. operators 
Computing M. For any v,w ^ T^T, let u = -u o u""^ and w = w o we have 

2 _ /" „V7„ V7„ „IV7„ |2 



^{u){v,w)v ■ W = IIuiv,v) ■ IIuiw,w) - IIuiv,w) = / pVpv,v^Pw,w - p\yPv,w\ dx. 

Assuming that the hypersurface St and v are sufficiently smooth we single out the leading order 
term of ^{u){v, ■)v. 

We first estimate IIu{v,v) ■ IIu{w,iI)). From the Divergence Theorem and (|3.11|) . 



JR"\St J St JR"\St 
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Using H3.6|) and (|3.1U|) for the above boundary integral and applying the Divergence Theorem twice 
for the interior integral, we obtain 

JR^-^St JR"\5t JSt + 

+ U_{wl,wl) + Vw^w+ ■ iV+ + Vyo.w- ■ iV-} - /9+w+Vtt,_^p+„ - p-W^Vw_p~^ dS 

Using (|3.13p to compute S/w±w± ■ N±, we obtain 

pVpv,v'^Pw,wdx\ < C(['u;||2(5^) + \w\hiM"^St)y 



1/ 



To compute IIu{v,w)'^, we use the decomposition = ■^'H±p^.^ — Aj}tTi:{DvDw) and the the 
Divergence Theorem to obtain 

[ p\Vp^,yj\'^dx = [ p^^Mp^^dS+ [ p\VA^hiiDvDw)\'^dx. 

iR"\5t JSt ' ' iR"\St 



Therefore, 



p\'^Pv,w\ dx- / Pv,w^Pv,wdS\ <C\w\^2(^n^St)- 

J St 



Moreover, from expression H3.14() of p^^^, we claim the terms other than A/" • (w^vj.) are of 
lower order. In fact, since 

and it is divergence free, its normal component on St is in H^^{St) and we have 

\piu,+^'^T^ ■ KK - ^~))\Hh(^St) - 
Therefore, summarizing these estimates, we obtain 

(3.24) mu){v,w)v-w+ [ \M--2V-{wi{vl-vI))\^dS\<C{\w\l2^St) + \<HR--.St))- 

J St 

Using Lemma I3.1L we define a self-adjoint positive semi- definite operator ^o{u){v) on TuT which 
is the leading order part of ^{u){v, ■)v. In Eulerian coordinates, ^o{v) takes the form 



(3.25) 



f± = ^n±M-^M^V,j_,-rN''V . {wi{vl - vl)). 



p+p- 

Assuming smooth St and from ()3.24|) . we can write 

^{v,w)w = iMq{v)w + at most 1st order diff. operators. 

Clearly ^o{v) is a second order negative semi-definite differential operator. Therefore, the linearized 
Euler's equation 1)3. 19(1 would be ill-posed if there had been no surface tension, for Mq{v) would 
become the leading order term. This is the Kelvin- Helmholtz instability of vortex- sheets. 

Note that since £^ is positive definite and is higher order than the self-adjoint ^q{v), it is not 
difficult to see that the linearized problem (|3.19p is well-posed. For a priori estimates of the Euler's 
equation (E) with the boundary condition (BC), the positive semi-definiteness of the leading order 
part £/{u) of &'^S{u) suggests to consider the inner product of (|3.5() with {Si'^S)^ut to obtain a 
priori estimates. 
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4. Main Results 

In this section, we will derive local energy estimate. We show that solutions of (E) with boundary 
condition (BC) are locally bounded in 

(4.1) v{t, •) G m''{M.'' \ St) and St G 

where k is an integer satisfying |A; > ^ + 1 (equivalently |A; > ^ + |). 

Definition of the energies and statements of the theorems. The conserved energy of the 
Euler's equation is given by Eq = f^n^g ^\v\'^dx + S{St). Higher order energies are based on the 
linearized Euler flow and thus involve the differential operators £/ defined in (|3.21j) and &t defined 
in dHSI). 

Let io^ : M" — > M", often written as to for short, represent the curl a vector field v defined on 
W \ St, i.e. 

uj{X) ■ Y = Vxv ■ Y - Vyv ■ X 
for any vector X,Y £ M". Viewing w as a matrix, its entries are ujj = uj{-^) ■ = diV^ — djV^. 

Definition 4.1. Given domains with i7+ compact and the interface S in H2^^^ and any vector 
field V G i?2^(M" \ S) with + V-\s = and V ■ v = 0, define the energy E{S,v), often written 
as E for short, 



I. A ,0 1, 



3 



where is the pressure due to the surface tension defined in (|3.17)) . 



Since the free boundary is evolving, we consider the following type of if 2 2 neighborhoods of 
hypersurfaces to maintain uniform constants in the energy inequalities . 



Definition 4.2. Let A = A(5, |fc — be the collection of all hypersurfaces S such that there 
exists a diffeomorphism F : S ^ S <Z , with \F — ids\^3^_i^^^ < S. 

Fix < 5 <C 1 and let Aq = A(5o, f A: - i, 6). From (IXTTl) . dTTrT) . and ^Ml, 



(4.2) b-lH»+^(M.^5) ^ ^\''\hHS), s G [i, h - i] 

3 3 

(4.3) KlL(//«(K"^5),//=-3{n)) <C, s G [4 - -k, -k - 1] 

where C is uniform in S G Aq. The next proposition gives bounds on the velocity and mean 
curvature in terms of the energy E. 

Proposition 4.1. For S G Aq with S G H^^^^ , we have 

for some integer m > depending only on k and n and some constant Cq > depending only on 
the set Aq. 

The proof of this proposition will be given below. Using this result we will prove the following 
theorem on energy estimates. 

Theorem 4.1. Fix 6 > sufficiently small. Then there exists L > such that, if a solution of (E) 
and (BC) is given by St with St G ^1''+^ and v{t, •) G C^{H^^{W'' \ St)), then there exists t* > 0, 
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depending only on \v{0, ')\fj^k^-^„ ^ y L> '^'^^ ^he set Aq, such that, for all t € [0, t*], 



(4.4) 



St E An and \k\ 3^. < L, 

E{St,v{t,-)) <2E{So,v{0,-)) +Ci+ [' P{Eo,E{St>,v{t',-)))dt' 

Jo 



where P{-) is a polynomial of positive coefficients determined only by the set Aq and Ci is an 
constant determined only by \v(0,-)\ 3,, 3 , and the set An. 

Since the domain is evolving, the above continuity assumption of f in t means that there exist 
extensions of both and V- to [0, T] x M" which are continuous in i^2^(M"). 

In order to prove Proposition 14. II and Theorem 14.11 we need the fohowing lemmas. 

Lemma 4.2. For any 5 G Aq with k G H'^{S), s G — |, |A; — 1], we have 

\^\h^{S) + I^Ih=+i(5) < C{1 + |k|h''(5))) 
for some C > uniform in S (z Aq. 

Proof. The proof is a straightforward application of elliptic regularity applied to 

(4.5) -AsU = -V^K + {\U\'^I - kU)U. 

and the fact that IIIL < C uniform in S" G An. □ 

Corollary. Suppose S £ Aq with k G H^''~'^{S), g G H^'^~^{Q^), and q = —A'^^g, then we have 

for some C > uniform in S (z Aq. 

The proof of this corollary follows Lemma 14.21 and the identities 

V7v±V7v«^g = M{N±) • Vq + D\{N±,N±) = M{N±) -Vq-g- k±V7v±9 on S 

- AVn^^ q = VNn^g-2D^q- DNn^ in 0± . 

Lemma 4.3. Suppose S" G Aq with n G H^^^^^S), 

1 13 3 1 

\{-As)'2 -M±\j,„s>,q^-,<C(l+\K\ 3,. 3 ), S' G [ k,-k 1. 

Proof. From the identity 

{-As-Nl)f = K±M±if)-2VNj-A±)-\DNn^ ■ fn^) - M±{N^) ■ (A4(/)iV± + V^/) 
for any smooth / : S" ^ M, and lemma IX2I we have 

I - As - ^hbiH"' {S),H'''-^(S)) - 

Using commutators estimates [— A5,7V±] and factorization, the lemma follows. A detailed proof of 
the commutators and the factorization is given in section 6 of |SZ06 . □ 

Proof of Proposition HTTl The two terms \'^^v\'j^2(^g-^ and |j2/*^^2 Vpk|^2(5) can be written ex- 
plicitly using the definition ()3.21() of ^ 

(4.6) \^-'v\hm^^s) = f vi{-AsM)''-H-^s)vidS 



s 



|2 

Il2(R"\5) 



(4.7) K^"5V^5«|i2«„^5^ = / K+M{-AsM)''-^K+dS, 
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where 

(4.8) M = {—M+)M-\—M-) = ((— AA+)-i + 

P+ P- P+ P- 

Clearly M is also self-adioint and positive. The estimates on IkI 3,. , and Ifil Si, 1 fohow 

immediately since from lemma U3lA/± behaves like (— A5)2. 

To bound v in terms of £"0 and E we note that Aw is controlled by E from 

(4.9) Aw* = djuj]. 

Therefore it is sufficient to control the boundary value V n± v± by E. Moreover from the identity 
Vjv± v± = [Dv±)*{N±) +ijj±{N±) where u}± is the restriction of uj on S", it suffices to show that Eq 
and E control v± = {Dv±Y{N±). 

We first estimate vj_ using the identity 

(4.10) A5 ul = V^{V- vl) + Ric((Z)r;+)*(iV+)T) + [Vx, cuJ^)(X,-), at :e G 5 

where Ric is the Ricci curvature of S, {Xi, . . . is any orthonormal frame of TxS, and the 



tangential curl uJ of uJ, is defined as 



a;T (x)(X) -Y = Vxi^I-Y-Vyi^I-X = Vxi^I-Y-Vyi^I-X 
for any X,Y G T^S. From the definition of i'± we have 

u;J^(x) ■Y = u+{x)- Vy v+ - n+(y) • Vx v+ 

Therefore, by Sobolev inequalities, there exists C > uniform in S* G Aq so that 

• I I 3fc 5 + I Ric I 3fc 5 l-^+l 3fc 3 IDV^I 3fc o + [tJ,, , I 3fc 5 
<|2?-l/X| 3 I. 5 +C|V+| 3,, 1 , . 

Here the norm H^^~a is chosen to illustrate that the term is lower order. In fact any H^^~°' with 
< a < ^ works. To estimate the divergence term V ■ uj^, one may compute 

V-ul = Asvi-V-{U+{vl)), 
which along with Lemma 14.21 implies 

YD ■ uL\ 3k 5 <Cbi| 3,, 1 +C|n+| 3,, 3, \v+ — (v+ ■ N+)N+\ 3,, 3, , 

<C|vi| 3,, 1 +C(1 + |k+| 3,, 3, )\v+\ 3,, , 

Therefore, from (|4.6|) and (|4.7|) . we obtain 

(4.11) ^ + + 

Finally, we only need to estimate = S/n+v+- A^+, much as in the way in jSZOGj . Extending 
into 17+ as = {Dv+)*{N-}-(^), where Nt-i^ is the harmonic extension of into Q"*", and comparing 
the two ways of computing V • on S" using 1) frames and uj and 2) divergence decomposition 
formula on S, we obtain 

1/+ • iV+ = iVx,u;)iX,) ■ N+ + Dv+ ■ DNn+ -V-uJ- k+u^. 

Moreover, 

VN+y+ ■ N+ =Vjv+(Vjv„^?;+ ■Nn+)-N+- Dv+{M+{Nn+)) 

=N+{vi) + V^+A;1A(Va.^^ v+ ■ Nh+) - N+ ■ Dv+{N+{Nn+)). 
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Therefore, from the estimate on D • i/J, we obtain 

\M+ui\ 3k 5 < CE^ +C(l + E^)\v:\ 3, 1 
which, along with ()4.9p and (|4.1ip . implies 

|2 



The estimate in Proposition 14. II follows immediately from Sobolev inequalities. □ 

Proof of Theorem 14.11 To prove Theorem 14.11 in addition to Proposition 14.11 we need the 
following: a) the estimates on the Lagrangian coordinates map and consequently k G H2^~2(^S), 
b) estimates on u; = Dv — (Dv)*, and c) commutators involving D^. In the following all constant 
C > will depend only on the set Aq. 

Estimate of the Lagrangian coordinate map u{t, y). We will only work on the domain Q\. From our 
assumption on v, the ODE ut{t,y) = v{t,u{t,y)) solving u is well-posed. Since u{t,-) : 
is volume preserving and |A; > ^ + 1, it is easy to derive, Therefore, 

(4.12) \u(t,-)-I\ 3^, .<C I \v{t\-)\ 3,, , lM(t',-)|^l dt\ 

where C > depends only on n and k. Let /i > be a positive large number specified later, 

(4.13) to = sup{t I \v{t\ Ol^i,^^,^^^^^ < /X, Vt' G [0,t]}, 

We have to > due to the continuity of v{t, •) in H'^^{nt)- From ODE estimates, there exists ti > 
and C2 > which depend only on /i such that, for all < t < min{to,ti}j 

(4.14) Kt,.)-/|^3,^^,^<C,t. 

It implies the mean curvature estimate, for all < t < min{to,ti}, 

(4.15) |.(t,.)|^3,_.^^^^<|.(0,.)|^3,_.^^^^^ + C73t. 

Here it is easy to see from local coordinates that C3 is determined only by /i and the set Aq. 
Therefore, there exists ^2 > determined only by ^ and the set Ao such that St G Ao for 
< t < min{to,t2}. 

Evolution of the curl uj = Dv — (Dv)* . It is easy to compute 

BtLO = DBtv - {DBtv)* + {{Dv)*f - {Dvf = {{Dvff - {Dvf = -{Dv)*uj - ujDv. 
Since |fc > S + 1, we have 



d 



(4.16) — / \D-2^-'voY 3^ , dx<C\v\ 3^, , 

The commutator involving D^. First we list a few commutators calculated in |SZ06j : 

(4.17) nt'^g = V-Dtg-{Dvr{Vg), 

(4.18) [Dt, A±i]5 = A^H'^Dv ■ D^A^^g + Av ■ VA^^g), 

for any smooth function g. It also has been proved in |SZ06j that, for any function / defined on St, 

7 3 3 1 

(4.19) \[Bt,^sMLiHnis^),Hn~HS,)) < ^I^L§*.K"..S.) e {- - ,,k,-k-- 



(4.20) l[Dt,A4]|L(//»2(so,H=2-i(SO) <C|^'U4fc™„, c^ ^2 G U,i^k--]. 



H^^iw^-^St) '2 2 2 2 

H5«(R"\S't) ^ '2' 2 2 
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It also implies that 



(4.21) 
(4.22) 



V\ „3, 



-St) 



.13, 3^ 

s G --, -k 

^ 2'2 2^ 

.1 3, 1, 
4'2 2^ 



Evolution of E. For the rest of this section, let Q = Q{\v 



K -3, 



positive polynomial in Sj. 



and \k\ 



denote a generic 



with coefficients depending only on the set 



Aq. One may notice here, by 

Lemma|1121 |n| s,, , , and If i I 3,, 1 , and IfTl Si. 1 can also be included in Q. Recall 

that 



Since 
then 



Bt^dS = {V-vl + K±v^)dS. 
■ v±\st = V ■ vl + K±v^ + Vn±v± ■ iV± = 



\k±v^ + V- v±| 3,_3 = \Vn±v± ■ N±\ 



<C\v\ 3,, 



and thus T)tdS would not complicate the estimates since |fc > ^ + |. 



Before we embark on calculating the energy inequality it is helpful to recall two facts. First, we 



only need to keep track of terms which can not be bounded by \v±\ 3^ and |av| 
second, Af± are selfadjoint operators of order one. 



1 , and 



2 dt I 



n+Mi-As,M)''-^'Dt+ K+dS 



St 



< Q. 



From dmH), and it is clear 

Id, 1 ^ ,9 

2dt' 

From the expression for Dt^K-t given in 

(4.23) Dt±/t± = -A5,4 - vi\Ii\^ + (P • n±)(4) = -As,fi - 4|n|2 + V„TK± 
and Lemma l4. 21 we only need 

(4.24) I / K^M{-AstMf~^^^TK+dS\ < Q. 

Jst + 

to derive estimate I. By considering a flow </>(r, •) on generated by Ti+v J, the above commutator 
estimates applied to D,- allow us to pull V^t to the front with lower order errors 



J St J St 



K-L I dS 



< Q. 



Therefore, inequality 1)4. 24() follows from the Divergence Theorem and I follows consequently. 



II: 



iill^^^lhiR^^St) - + Is, vi{-As,M'fK+ dS 
where the extra term 

P+ 



< Q 



E,. 



2(p+ + p_) Js, 



V X K+ ■Mi-As.MT-'V^-T K+ dS 



P- 



2(p+ + p_) Js, 



\7,T^+-Mi-As,M)'"''V,TK+dS} 
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From (gSl), (EHH), and it is clear 



(4.25) 



^ d , .k 
- — \M2v\ 

2dr 



St 



< Q. 



Using (prra . dsn), (EHI), and (jSHJ), we have 

^ N+pI, - VN+A-hr{Dvf-MK+ - V^Tvi + U+{vl,vl). 



P+ 

From the corohary of Lemma 14.21 we have 

(4.26) \VN±A±hT{Dvf\ 

Therefore, 



(4.27) 



1 1 I 
' 2 (it ' 



< Q. 



1 



St 



v'^{-l:^StW-\-^St){ - — AT+p^,. - A/'K^ 

P+ 



-y^Tvi + Ii+{vl,vl))dS\<Q. 
From ()3.1ip . we can write the first term in the second integral as 

-—M+pI, = —M+N-\2V^-r Tvi-Ii+{vl,vl) - Ii^{viy_) 
P+ P+ + " 

- V N+A+^ti{Dvf - V7v_AlHr(Dw)2} 

Since \M^^J\f-L.\ 3,1 < C, we have that its dual AZ+AA^^ satisfies the same estimate. There- 
fore, using (|4.2H|) . we obtain 

I - —M+pI, - —M+M-\2V.Tvi-Il+{vl,vl)-U^{viy_))\ , . < Q. 

From Lemma l4.,'S| we have 



-1 P+P- 



p+ + p_ 'L{Hi'=-i(S't),H§*-i(5t)) 



which, along with Lemma l4.2( implies 



1 



—N+p- 



P+ ""^ P+ + P 
Substituting this estimate into 1)4. 27() . we have 



(2 V,T _,T vi-Ii+{vlvl)-Ii^{vL,vL))\z 



H^"'-^(St) 



< Q- 



\- — \ 



2V\ 



v^{-As,M)'-\-As,){ 



St 



P+ + P- 



U+{vl,vl) 



P- 



P++P- 



2p_ 



^P+ + P- P+ + P- 
Much as the proof of inequality ()4.24|) . by commuting V^t we obtain 

2p. 



vI)}dS\ < Q. 



St 



Therefore, we have 



vi{-As,M)'-\-As,)Vvi ■ {P—^v. 



P+ + P~ 



P+ + P^ 



-vl)dS\ < Q. 



(4.28) 



Id, k 2 
\2'^\'^^^\L2{R"-^St) 



[ vi{-As,M)''-\-As,){ 

J St 



p+ 



P+ + p- 



p- 



P+ + P- 



-U^{vl,vl) -MK+}dS\ < Q. 
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In order to estimate the terms with 11, we first use Lemma 14.21 and identity ()4.5|) to obtain 

\As,iU±{vl,vl)) - 'D\±{vl,vl)\3,_s^< Q. 



Since 
we have 



V\±ivl,vl) = V,tV,tk± - V^Tvl ■ Vk± 



+ ^~ \7^TV^TK-}dS\ < Q. 
P+ + P- 

As in the proof of H4.24|) , we apply the commutator estimates (|4.19l) and (|4.2U|1 to move one of V„t 
and obtain 

P+ + P- J St ~ ~ J St 

From identity 

which imphes 



/ V,T-Dt^K+-M{-AstAr)''-^V,TK+dS+ [ vi{-AstM)''K+dS\<Q. 

P+ + P- J St ~ ~ J St 

Using ()4.19() and H4.2U() one more time, we obtain inequahty II. 



Proof of Theorem \4.1\ Adding inequahties 1)4. 16() . I, and II, we have 

(4.29) E{t) - E{0) + Ee^t) - E,.M < ^ Q{\v\HZ^(^^n^s,). I^Ih3'=-i(s,,)) 

where Q is a polynomial with positive coefficients that depends only on Aq. This inequality holds 
on [0, min{to, ^2}] where to is defined in (|4.13|) and t2 is determined only by |t>(0, ■)lH3'=(M"\St) ^'^^ 
the set Aq, the neighborhood of Sq in H2^~2. Clearly 

I eXl - I |^3fc_b^j^^^^^^| |^3^_3^^^^, 

where C > depends only on the set Aq. Interpolating v between if^^ (M" \ St) and H^''~^ (M" \ St) 
and K between H^''~^{St) and H^''~2(^St), we obtain from Proposition 14.11 

1 



\Eex\ < -E + Ci{l + \v 



3 



2 '//•^''-^ (R"\5t) 

for some integer m > where the constant Ci, which include | determined only by 

H ^ (St ) 

£"0 and the set Aq. From the Euler's equation (|;-{.18|) . (|4.2j) . LemmaE21 and its corollary, 

We can use the Lagrangian coordinate map u{t,-) to estimate 
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Through a similar procedure of the derivation of H4.14() . there exists > 0, depending only on 
|f(0, •)l_H'3fc(iKn^5j-) and the set Aq so that for < t < minjiojis}, 

< [ Qdt' 
Jo 

for some polynomial Q with positive coefficients. Therefore, 

where Ci is determined only by |f(0, •)! 3 and the set An. Thus 

E{St,v{t,-)) <2E{So,v{0,-)) + Ci+ [ Qdt'. 

Jo 

By inserting the above inequality into (|4.29|) and using proposition 14. H we obtain (|4.4j) . By choos- 
ing /i large enough compared to the initial data, Theorem 14.11 follows . □ 

Notation 

A*: the adjoint operator of an operator. 

D and d: differentiation with respect to spatial variables. 

V/: the gradient vector of a scalar function /. 

Vx- the directional directive in the direction X. 

_L and T: the normal and the tangential components of the relevant quantities. 

= (9t + v'^dr^i the material derivative along the particle path. 
St = dQt the boundary of a smooth domain evolving in time. 
N{t,x): the outward unit normal vector of St at x £ St- 
11: the second fundamental form of St, Il{t,x){w) = £ TxSt- 

U{X,Y) = U{X) ■ Y. 
k: the mean curvature of St, i.e. k = trll. 
/t^ = the harmonic extension of / on ^If 

J^{f) = ^ N'H{f ) : 5 ^ M: the Dirichlet-Neumann operator. 
X = X o u"^ the Lagrangian coordinates description of X. 
D: the covariant differentiation on St C M"". 
T^w = for any x e St we T^St- 
TZ{X,Y), X,Y G TxSt- the curvature tensor of St- 

— trP^: the Beltrami-Lapalace operator on a Riemannian manifold Ad. 
A~^: the inverse Laplacian with zero Dirichlet data. 

r = {$± -.nf ; volume preserving homeomorphism, such that = <I>_(9$7 )} 

^: the covariant derivative on T, 

Si: represent S in Eulerian coordinates. 

the curvature operator on F. 

represent M in Eulerian coordinates. 
II: the second fundamental form of T C 
II„('u;i,u;2) = ^w^W2, for any u G T, wi,W2 G TyT 
Pv,w = —A~^tic{DvDw). 
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